Interferometry using spatial adiabatic passage in quantum dot networks 
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We show that techniques of spatial adiabatic passage can be used to realise an electron interferom- 
eter in a geometry analogous to a conventional Aharonov-Bohm ring, with transport of the particle 
through the device modulated using coherent transport adiabatic passage. This device shows an in- 
teresting interplay between the adiabatic and non-adiabatic behaviour of the system. The transition 
between non-adiabatic and adiabatic behaviour may be tuned via system parameters and the total 
time over which the protocol is enacted. Interference effects in the final state populations analogous 
to the electrostatic Aharonov-Bohm effect are observed. 

PACS numbers: 05.60.Gg 73.63.Kv 73.23.Hk 



I. INTRODUCTION 

Quantum information science offers a wealth of new 
and potentially important phenomena to explore. One 
hope is that these possibilities will translate into prac- 
tical and unique devices, with functionalities unmatched 
by classical devices. One field that has explored quantum 
mechanical effects and quantum coherence in particular 
is the field of light-matter interactions, and especially 
the couplings and effects investigated in the physics of 
optically driven multi-level atomsi. Moving from atomic 
systems, with particular level structures, to spatial sys- 
tems, perhaps defined lithographically, provides dramatic 
opportunities to tailor the Hilbcrt space and connectiv- 
ity of the resultant quantum system. This allows natural 
translations of known effects from the atomic realm, to 
the electronic realm. Examples include charge qubits 2 
and qudits 3 -, adiabatic passage^, and electromagnetically- 
induced transparency^. A review of some of the progress 
in mesoscopic electronic systems, focussing on the quan- 
tum electronics/quantum optics/quantum information 
interface can be found in Ref£. 

The solid-state provides exceptional possibilities for re- 
alizing spatially defined Hilbcrt spaces, especially in the 
context of superconductors, donors in silicon! and quan- 
tum dots. Here we focus on the latter for conceptual 
simplicity, but our ideas are equally relevant to all these 
implementations. Quantum dots provide a rich platform 
for exploring novel experiments in quantum mechanics 
and often allow for fine control over many system param- 
eters. Chains of connected dots are of particular inter- 
est, with recent studies including examinations of triple 
quantum dots in GaAs 2DEG structures^ 
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car- 
bon nanotubesi 3 - and double quantum dots in silicon^, 
which extend the accessible spatial Hilbert space to allow 
clear connections with the quantum optical case. 

Coherent Tunnelling Adiabatic Passage (CTAP)^ has 
been proposed for transporting quantum information and 
is a spatial analogue of the well known STIRAP protocol 
from quantum optics^. It transports a particle coher- 



ently using a counter-intuitive coupling sequence of tun- 
nelling matrix elements. CTAP has been recently demon- 
strated by Longhi et al. using photons in three and multi- 
waveguides structur o 15 ' 16 ! 17 but it has also been proposed 
for observation in quantum dot a 18 ' 19 i 20 . Other platforms 
for which demonstrations of CTAP have been proposed 
include phosphorus donors in silico n 21 i 22 i 23 , where CTAP 
can be used as the transport mechanism in a quantum 
computer architecture, and also in superconductors 2 ^, 
single atoms in optical potential a 25 i 26 , and Bose-Einstein 
Condensate a 27 i 28 . STIRAP has also been discussed for 
spatial particle motion in quantum dots, see for example 
Ref,« 

The extensible nature of controlling the spatial loca- 
tion of states and their connectivity has led to CTAP 
being extended to multiple recipients (MRAP) 31 i 32 as 
a means of implementing a form of fanout for a quan- 
tum computer. This branching ability, and the fact that 
CTAP is a coherent process, lends itself to investigat- 
ing interferometry with devices using the CTAP as the 
transport mechanism. 

Interferometry is a well-known means of probing the 
wave- like properties of particles, and non-trivial quan- 
tum phases. There have also been solid-state experiments 
that utilise interferometry to observe the wave- like nature 
of electrons in Aharonov-Bohm (AB) and related inter- 
ferometers for 'which-path' measurement a 33 i 34 ' 35 . The 
AB effect is a topological effect where a charged particle 
(for example an electron) traverses a loop. If the loop 
encloses a magnetic flux, this will break the symmetry of 
the paths and introduce a non-trivial phase. Interferom- 
etry then allows this phase to be revealed in oscillations 
in the final state population, as a function of the enclosed 
flux. 

Here we are considering an adiabatic realization of the 
electrostatic Aharonov-Bohm (EAB) effect^^. This ef- 
fect is named by analogy with the AB effect, but is not 
a topological effect. Instead, an external electrostatic 
field is used to break the symmetry of the two paths, 
which again gives rise to non-trivial phases and popu- 
lation oscillations as a function of the field. The EAB 
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FIG. 1: The interferometer structure consisting of six coupled 
sites with gate-controlled TMEs Qi(t) and Q,2(t). The particle 
is moved from |1) to |5) using the alternating CTAP protocol 
via the two intermediate states |3 U ) and |3d). The control 
gates, labelled A u and Ad, break the symmetry of the paths, 
leading to interference effects seen in the final state population 
as a function of the magnitude of the energy shifts and the 
total time of the transport. Note that we assume that the 
energy that the electron has at each site has been set equal. 



device we consider uses adiabatic-passage as the means 
of transport. This provides interesting flexibility and a 
rich phase space to explore. In particular this six-dot 
geometry is the simplest topology that permits both adi- 
abatic passage and a non-trivial loop. Furthermore, our 
model shows an interesting interplay between adiabatic 
and nonadiabatic features, as is discussed below. 

In our proposed quantum dot geometry, the Tun- 
nelling Matrix Elements (TMEs) are controlled via sur- 
face gates using the alternating coupling scheme with 
counter- intuitive pulse ordering (described below) , a vari- 
ation of the CTAP protocol called Alternating CTAP 
(ACTAP). This geometry has been investigated in linear 
chains for quantum dots^ and in the context of a five 
donor Si:P device 4 ^. Again, we do not wish to explore 
the microscopic details of particular implementations too 
strongly, instead focussing on the physics that can be 
demonstrated in all potential platforms. Because of the 
wide range of possible implementations, we do not treat 
decoherence here, although its effect will be to wash out 
the interference patterns. 



effect has been investigated in the context of one dimen- 
sional mesoscopic rings in metallic* 3 ^ and semiconductor 
structure a 39 i 40 and generalised for the case where both 
electrostatic and magnetostatic potentials exists. 

We show that CTAP can be used to explore physics 
similar to those seen in the EAB effect. Our suggested 
geometry, highly reminiscent of traditional AB rings, is 
shown in Fig. [1] It consists of six sites with the tunnelling 
between sites controlled via surface control gates. The 
energies of the two middle sites in each branch may also 
be controlled by a gate and it is assumed that the energy 
of an electron at each other site is also controlled, static, 
and equal. A square quantum dot, such as that examined 
in Refsi 42 i 43 , with added coupling to a source and drain, is 
one possible implementation of this structure. Aharonov- 
Bohm-like oscillations in an adiabatic passage four dot 
ring in the presence of a magnetic field has also been 
considered 4 ^. 

The device shown in Fig. Q] is analogous to a Mach- 
Zehnder interferometer where the final state population 
observed is dependent on the interference condition cre- 
ated by the gate-controlled energy shifts on either arm. 
Because of the deep analogy between this structure and 
that of quantum optical systems, we shall refer to the 
energy shifts of these donors as dctunings. The detuning 
on the middle sites is controlled individually by gates and 
plays a similar role to the phase difference between the 
two arms of the Mach-Zehnder interferometer. Here we 
are investigating the passage of a single electron travel- 
ling through the device in some superposition of the two 
pathways. Where the detuning is different in the two 
paths there will be different phase accumulated by the 
electron, which is exhibited in the interference pattern in 
the final state population of the end-of-chain site. The 



II. ANALYTICS AND MODELLING 

To explore the dynamics of the ACTAP interferometer, 
we write the Hamiltonian for this system as: 

H = n 1 (|l)(2| + |3„)(4| + |3 (i )(4|) + / 1 .c. 
f2 2 (|2)(3„| + |2)(3 d | + |4)(5|) + /i.c. 
+A„|3„)(3„|+A d |3 d )(3 d | (1) 

where f2i and ^2 are the gate controlled tunnelling ma- 
trix elements (TMEs) and A u , Ad the energy detunings 
of sites |3„) and |3<j) respectively and % has been set 
to 1. In general the TMEs between nearest neighbours 
will be independent, however by construction, we have 
chosen gate values that ensure the form of the Hamilto- 
nian shown in Eq. 1. This form is necessary to enforce 
the ideal ACTAP coupling scheme. It should be real- 
ized, however, that complete symmetry is not required. 
The ACTAP scheme has some robustness to variations 
from the ideal case, and this was explored for the linear 
ACTAP chain in Ref4£. The choice of equality greatly 
simplifies the analysis and encapsulates all of the physics 
of CTAP. Indeed this is the simplest coupling scheme 
that provides intcrferomctry with a non-trivial loop and 
CTAP coupling. For simplicity, we choose squared sinu- 
soid functions for the form of the TMEs: 



^1 = ill max sin 2 ( — ) , 

^2 = ^2 max COS 2 ( — ^— ) , (2) 

where Qi max = ^2 max = ^ max is depicted in Fig. [^a). 
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FIG. 2: (a) Tunnelling matrix elements Qi(t) and fla (t), 
effecting the counter-intuitive pulse sequence, (b) Eigen- 
spectra with varying detuning, A„ = Ad = (solid lines), 
A» = -A d = 0.25n max (dashed) and A u = -A d = Q max 
(dotted). The eigenstates can be ordered as indicated. The 
double degeneracy of the null space is lifted with nonzero 
detunings of the central sites, we arbitrarily label the more 
energetic eigenstate of the pair |23q + '), and the less energetic 



In the case A u = Ad = the time dependent eigenval- 
ues of the Hamiltonian are 



E 
E± 

E 2 ± 



= 0, 



± 



30 2 + 302. - yjnf + wofo 2 . + o| 
71 ' 



'30? + 302 + ym + iir2?oTTo! 

: — — — ' (3) 



where E is doubly degenerate. These eigenvalues are 
plotted in Fig. [3Jb), along with eigenspectra with non- 
zero detuning for selected values of A„ = — A<j, for com- 
parison. The corresponding states are labelled \Vq), 

[Dq ), |2?±)i an d l^2±)- They are in general difficult 
to represent in closed form, however the eigenvalues or- 
der naturally, and these are indicated in Fig. HJb). In 
certain limits we may extract some useful information 
about these states, and some of these will be given be- 
low. These plots illustrate the lifting of the degeneracy 
of the \T>q) states with detuning, in these cases with op- 
posite detunings, ie A u = — Ad- 

To enact the CTAP protocol we initialise the device so 
that the particle occupies site |1) at t = 0, and apply an 
alternating pulse sequence in the counter-intuitive order- 



ing along the chain to effect population transfer. Con- 
ventional CTAP in linear devices (e.g. the straddling 4 ^ 
or alternating scheme o 18 i 45 ) works by adiabatically trans- 
forming the null state from this initial state to the desired 
output state. Although in general there is no null state in 
our scheme, the adiabatic transfer still works analogously. 
The evolution of this state under the counter-intuitive 
pulse sequence gives rise to a smooth change in the pop- 
ulation of the sites from |1) at t = to |5) at t = t max . 
Since we are using the alternating coupling scheme there 
is a transient occupation of |3 tt ) and \3d) during the pro- 
tocol. In the straddling scheme 4 ^ occupation of all the 
intermediate sites along the chain is strongly suppressed 
limiting its effectiveness for interferometric sensing. Note 
that as the degeneracy of the central sites, |3 U ) and |3<j), 
is broken by the gate induced detunings, interference will 
be observed in certain regions of phase space. This in- 
terference is the effect that we are seeking to understand 
and exploit in this work. 

To examine the effect of the detunings on the final 
state population for the six site system we numerically 
solve the Schrodingcr equation in Eq. (TjTJ for the final 
state population of |5) after the CTAP pulse sequence 
has been performed in a finite time, long enough to en- 
sure adiabatic transfer in the zero detuning case, as a 
function of the detunings on the middle sites. The re- 
sults of this calculation are shown in Fig. [3] This map 
of the final state population (p^s) shows several regimes 
of interest. There is a prominent 'cross' of high-fidelity 
transport where either A u = or Ad = 0. This region 
is governed by adiabatic timescales, and is discussed in 
Section III. Superimposed on this, is a line of alternating 
high- and low-fidelity regions where A„ = — Ad- This 
line corresponds to the EAB effect and is governed by 
nonadiabatic oscillations, and is discussed in Section [TV] 



III. ONE ARM WITH LARGE DETUNING, 
OTHER QUASI RESONANT 

The main adiabatic feature of the system mapped out 
in Fig. [3] is the central cross region of high fidelity trans- 
fer. It is most clearly observed where at least one of the 
middle site detunings is zero we have an adiabatic path- 
way that is similar to that of the simpler case with only 
one path as discussed i n 18 i 45 . To examine the width of 
the high fidelity region we consider the case where one 
detuning is taken to infinity as in this case we can re- 
trieve some relatively simple analytical expressions. In 
this limit we recover a five site system with detuning on 
the middle site, and we explore the dependence of the 
adiabaticity on the detuning of the central site. Since 
the system is symmetric, the width of the central cross 
feature is the same for the vertical and horizontal axes. 
We are not aware of this case being treated previously in 
studies of alternating coupling schemes for adiabatic pas- 
sage. It is more conventional to consider detunings for 
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FIG. 3: (Colour online) Map of the final state population 
(pss) for finite time showing the regions of adiabatic and non- 
adiabatic evolution. Dark regions indicate high fidelity trans- 
fer, the result of adiabatic evolution. Light regions indicate 
where the transport is low fidelity transfer. Note the overall 
'cross' of high-fidelity adiabatic transfer, with the EAB like 
oscillations along the line A u = — A^. 



states without any population (in our case the even num- 
bered sites), or the end of chains, and so this case rep- 
resents an interesting alternative detuning regime, which 
appears well suited to interferometric applications. 



Without detuning the eigenstates are eigenvectors of 
the five site chain can be fairly easily represented, as is 
done in Ref.— , but the central site detuning makes this 
problematic. However, we can write down the states in 
appropriate limits, namely at the extrema of the proto- 
col, and in the middle (when the adiabaticity is highest) 
as a series in the detuning, A. Note that for consistency, 
we shall keep the notation \T>q) for the CTAP transport 
state, and the other eigenstates similarly labelled, how- 
ever as mentioned above, there is no longer a null state 
with A ^ 0. 



At the beginning of the protocol, t = 0, we have Sli = and the eigenvalues and eigenvectors are: 

E = 0, E ± = ±Q 2 , E 2± = ±Q 2 -j + 0\A] 2 , (4) 

l»o) = |1>, (5) 

P±) = ^(|4)±|5», (6) 



^-^ +0[A] ^ |2) + |3) , CO 

/2±A +0 [A]2 

(8) 



and at the end of the protocol, t = i max , we have 2 = and 



E = 0, E ± = ±f2 l5 E 2± = ±n t - y + Q[A} 2 , (9) 
\Vo) = |5), (10) 
P±) = ^(|1>±|2», (11) 



±1-2H7+W J|3> + |4> 

|2>2±) = ^ 1 ' • (12) 

v /2^fT0[Ap 

(13) 

The more interesting case is at the midpoint of the transport protocol, i.e. t = t max /2, and for simplicity, setting 
Oi = fl 2 = ^max/2, we have 

E = ^ + O[A} 2 , £ ± =±^p, E 2± =±^p^ + ^ + 0{Af, (14) 

(|l)-|3) + |5))+(^ + 0[A] 2 )(|2) + |4)) 
\V ) = ^= 1 , (15) 
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(16) 

(17) 
(18) 



To quantify the adiabatic transfer in this five-site limit, 
and thereby gain insight over the adiabatic cross, we use 
the adiabaticity parameter defined for the five-site con- 
figuration 



A 



\E 4 



Eo\ 



(19) 



For adiabatic evolution of the system we require A <C 1. 
Applying a detuning to the central site the adiabatic- 
ity changes and may shift the system out of the adia- 
batic regime. Numerically one sees the increase in the 
adiabaticity parameter with increasing A as is shown in 
Fig. 0J High fidelity transfer (low A) occurs over a range 
of detunings. The effect of the detuning has been ex- 
amined for STIRAP- and for straddling CTAP in RefA 
However it is worth noting that these cases are different 
from the one we are considering at present, as they other 
treatments examine energy shifts of states with vanish- 
ing occupation. Here, the populations are explicitly non- 
zero due to the construction of the alternating coupling 
scheme, increasing the effect of the detuning. This is 
an important feature and in fact necessary to achieve an 
interferometric readout signature. With the sinusoidal 
pulse scheme, the maximal value of the adiabaticity oc- 
curs at t — ^1 max — ^2 ma: 

we finest 



A = 



-in 



\/3^max^max 



(20) 



To include the effect of the detuning on the adiabaticity 
we perform a series expansion on the adiabaticity in A. 
The adiabaticity parameter, to second order in A is 



A 



4tt 



20ttA 



V30 



max •rnax 



3\/3^ m ax^max 

56?rA 2 



(21) 



9 \/3^ ma x^max 

This form of the adiabaticity parameter is plotted in 
Fig. UJ and compared with a full numerical calculation 
of A. 



IV. INTERFEROMETRY IN THE SIX SITE 
SYSTEM 



Returning to the map of final state population as a 

function of the middle site detunings in Fig. [3J as well 

18 ~ 




FIG. 4: Maximum adiabaticity through the protocol as a 
function of middle site detuning for A-CTAP in a 5 site sys- 
tem, corresponding to the case in the interferometer where one 
detuning has been taken to oo. The solid line is the numer- 
ical determination, Eq. (|19|l and the dashed is the analytical 
result to second order in A, Eq. (|22p . The smooth increase 
in adiabaticity correlates to the smooth reduction in fidelity 
seen in Fig. [3] in the zones where A u and Ad have the same 
sign. 



as the central cross feature there is also the line of low- 
fidelity regions across the A u = — Ad diagonal where 
we see non-adiabatic evolution. Along this line we see 
regions of low fidelity transfer superimposed upon the 
region where we expect high fidelity. The position and 
number of these areas vary with the total time, as is 
shown in Fig. We see that as the total time is in- 
creased, the frequency of the nonadiabatic oscillations 
increases with increasing detuning. 



We can understand the interference by considering the eigenspectrum with non-zero detuning, in particular focusing 
on the lifting of the null state degeneracy. As we see in Fig. [2fb) when the sites are oppositely detuned (A„ = — A^), 

the states |I>q + ' > ) = \T>{ ') are degenerate at the start and end, but are split evenly during the protocol and it is the 
population oscillations between these states which will be important, rather than the adiabaticity treated in Eq. (1221) . 
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To first order in A, at the midpoint of the protocol we find 

4* 

The interference observed in the final state population 
can be understood as resulting from the electron under- 
going Landau-Zencr oscillations between these two states 
as the Hamiltonian is evolving with the protocol. The 
non-adiabatic behaviour can be thought of as arising due 
to the total protocol time required for the system to be 
able to resolve these two states. That is, when the total 
time is large compared to the energy gap between them 
the system is able to resolve the lifted degeneracy of the 
null states. This energy gap depends on the detunings 
of the central sites. Empirically, we find for a given time 
f max , maxima in the transfer fidelity occur when 



A„ = ^- (24) 



where A„ is the detuning of the n th maximum with 
A„ = ±A„ = T-A, and a factor / ~ 20 which has been 
determined empirically. Since the states we are inter- 
ested in are approximately parallel, a doubling in the 
total time will correspond approximately to a doubling 
in the frequency. Since these oscillations arc periodic the 
non-adiabatic behaviour will also be periodic in the total 
time. This is plotted in Fig. [5] where we see very clear 
linear relations on a log-log scale, except in the large A 
limit. With increasing detuning there is an overall re- 
duction in fidelity, as the barrier between sites 1 and 5 
increases, thereby suppressing population transfer. This 
effect is most evident in the lower fidelity transport re- 
gions outside the central cross feature in Fig. [3J 

The distinct regions of adiabatic and non-adiabatic be- 
haviour suggest that such a device might be interesting 
for charge sensing. Mapping out the rate of change of the 
final state population with respect to one detuning, as 
shown in Fig. [51 we can identify device detuning configu- 
rations which will display a large response to any change 
in the energy level of the one of the middle sites due to a 
charge being in close proximity. Using one arm as the de- 
tector arm and selecting a detuning configuration which 
will show this large response we can see that the final 
state populations respond to the presence of a near-by 
charge. Fig. [6](b) also shows the dependence of the sen- 
sitivity on time for the first fringe from Fig. [Ha) . The 
sensitivity of the interferometer at these operating points 
scales linearly with the total protocol time. 



(22) 
(23) 

I 




FIG. 5: (Colour online) (a) Final state population as a func- 
tion of total time for anti-symmetric detunings (A„ = —Ad). 
As the total time increases we may resolve more instances of 
interference when the accumulated phase difference between 
the two paths is it/2, (b) Trace of the final state population 
as a function of A u for the maximum time shown in (a). 



V. CONCLUSION 

Mapping of the final state population of a 6 site CTAP 
interferometer in an electrostatic Aharonov-Bohm style 
geometry shows an interesting interplay between adia- 
batic and non-adiabatic transport of an electron in the 
device. We may tune the device between these two 
regimes easily by modifying the total protocol time for 
the CTAP transport or by changing the on-site energy 
of the middle sites in the two different paths. We have 
modelled the behaviour of this device and the dependence 
of the detuning on the adiabaticity. We see effects simi- 
lar to electrostatic Aharonov-Bohm oscillations as seen in 
the periodic behaviour of the final-state population when 
the two paths have opposite detuning. In some config- 
urations a small alteration to the energy of a site will 
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FIG. 6: (Colour Online) (a)Derivative of the final state popu- 
lation with respect to A„ . Regions with a large rate of change 
may be considered useful for charge sensing applications. The 
very light regions are of high negative derivative and dark re- 
gions high positive derivative. The presence of a charge will 
alter the final state population due to its effect on the energy 
of the detuned site |3«) or |3<j) in the arm being used as the 
sensor, (b) Dependence of the sensitivity (derivative of the 
final state population) on i ma x. 



result in the device being shifted from adiabatic transfer 
to non-adiabatic. This is seen in a large change in the 
final state population. Such sensitivity suggests a device 
which may also offer opportunities for charge sensing ap- 
plications. 
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